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Abstract: We show that treating the blackbody radiation field as a heat bath enables one to utilize powerful techniques 
from the realm of stochastic physics (such as the fluctuation-dissipation theorem and the related radiation damping) 
in order to treat problems that could not be treated rigorously by conventional methods. We illustrate our remarks 
by discussing specifically the effect of temperature on atomic spectral lines, and the solution to the problem of 
runaway solutions in the equation of motion of a radiating electron. We also present brief discussions relating to 
anomalous diffusion and wave packet spreading in a radiation field and the influence of quantum effects on the laws of 
thermodynamic s . 
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1. Introduction 

Stochastic physics deals with fluctuations, principally ther- 
mal and quantum. The subject is often loosely referred to as 
"Brownian motion" since it was first seriously studied when 
Brown, in 1828, observed the random motion of pollen grains 
immersed in a fluid [1] at temperature T, in the absence of ex- 
ternal forces. Einstein [2] used a diffusion equation and showed 
that, for large times, the mean-square displacement is propor- 
tional to T/7, where T is the temperature and j^ 1 is the colli- 
sion rate, this being the first example of a fluctuation - dissipa- 
tion (FD) theorem. Soon after, Langevin [3] presented a simple 
phenomenological approach, by writing down the first example 
of a stochastic differential equation. In this equation, the total 
force acting on a particle due to its environment is separated 
into two parts: a frictional force and a fluctuation (random) 
force. These terms are very different in nature: The fluctuation 
term is basically microscopic in nature and has a time scale de- 
termined by the mean time between collisions whereas the time 
scale of the frictional force is proportional to the self-diffusion 
constant and is much larger. 

Later, another example of a fluctuation-dissipation theorem, 
arose in the analysis of so-called Johnson-Nyquist noise [4, 5]. 
All of this work was classical in nature but in 1951, Callen and 
Welton [6] presented a general quantum FD theorem. Such a 
theorem is implicit in the pioneering work of Ford, Kac and 
Mazur [7] who presented a microscopic quantum Langevin ap- 
proach to the case of an oscillator interacting with a heat bath 
composed of an infinite number of coupled harmonic oscilla- 
tors. This work was later generalized by Ford and Kac [8] and 
by Ford et al. [9] . In the latter paper, the earlier work was gener- 
alized by writing down what we referred to as the IO (indepen- 
dent oscillator) model, describing the system of a quantum os- 
cillator in an arbitrary potential, at an arbitrary temperature T, 
interacting with a heat bath of oscillators which were not inter- 
acting with each other. By assigning arbitrary masses and fre- 
quencies to all the oscillators, we obtained in essence a model 
which incorporated a variety of existing models. In particular, 
by means of a series of unitary transformations, we showed 
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that the blackbody radiation field (BBR) could be treated as 
a special case of the IO model. Concomitantly, because the 
Hamiltonian for an oscillator in a BBR field is universally ac- 
cepted, it provided in essence a "rosetta stone" in validating 
our choice of the more general IO model rather than any of the 
various models that one finds in the quantum optics literature. 
In particular, we showed [9] that the well-known RWA (rotat- 
ing - wave approximation) model of a heat bath has a serious 
problem in that the corresponding Hamiltonian does not have 
a lower bound. 

In this paper, we confine our attention, for the most part, to 
the BBR bath. Thus, in Sec. 2, we discuss the fundamentals un- 
derlying this subject and how they may be applied to the exper- 
imental results on the effect of temperature on spectral lines. 
In Sec. 3, we show how the quantum Langevin equation devel- 
oped in Sec. 2 forms the basis of a derivation of a solution to 
the problem of runaway solutions. In Sec. 4, we consider some 
miscellaneous applications. In particular, we show how mo- 
tion in a BBR field gives rise to anomalous diffusion and wave 
packet spreading, phenomena that are not apparently amenable 
to conventional QED techniques. Also, we point out that ther- 
modynamic concepts developed in Sec. 2 form the basis of ar- 
guments against the various claims that quantum affects could 
lead to violations of the fundamental laws of thermodynamics. 
In Sec. 5, we present our conclusions. 

2. QED shifts due to blackbody radiation 

Hollberg and Hall [10], using high-precision laser spectroscopy, 
measured photon heated Rb atoms to temperatures T as high 
1000-K" and analyzed the photon spectra associated with tran- 
sitions from the high Rydberg 36s state to the tightly bound 
5s state. They found an increase in photon energies propor- 
tional to T 2 which they concluded represented energy shifts 
due to temperature. Our conclusion [11, 12, 13] is that they 
have measured free energy shifts, as we will now argue. 

There is general agreement that the main frequency shift 
arises from T effects on the high Rydberg state since the ef- 
fect of T on the tightly bound state is negligible. In essence, 
we are dealing with a temperature dependent Lamb shift. Con- 
ventional atomic approaches to the problem have been carried 
out [14, 15] leading to the conclusion that the dominant energy 
shift ~ T 2 . The essence of the conventional calculation can be 
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simply obtained by an extension of Welton's T = calcula- 
tion [16] for the Lamb shift, and it also serves as a transparent 
foil to our approach. In this approach the weakly bound Ry- 
dberg electron is treated as a free electron which undergoes 
rapid oscillations due to the electric field associated with the 
BBR. The energy of oscillation W(u) of an electron moving 
in one dimension in an electric field Eoe~ lut is 



quantities for the heat-bath oscillators. Also f(t) is a c-number 
external force. The infinity of choices for the mj and u)j give 
this model its great generality. 

Use of the Heisenberg equations of motion leads to the QLE 
[9,11] 



W{uj) = e 2 E 2 /4muj 2 = (2ire 2 /3mu> 2 )(3E 2 /8ir). 



[1] 



+ / dt'fi(t - t')x(t') + V'(x) = F(t) + f(t), 



[5] 



Identifying 3i? 2 /87r with w(w,T), the energy energy density 
of the electromagnetic field, one substitutes the Planck distri- 
bution 



u{u,T) = (ftcj77rV)/[exp(?kj/fcT) - 1] 



[2] 



and integrates over all frequencies to obtain the mean energy. 
In three dimensions this is to be multiplied by a factor of three 
to give 



where V'(x) = dV(x)/dx is the negative of the time-independent 
external force and ^i(t) is the so-called memory function. F(t) 
is the random (fluctuation or noise) operator force with mean 

(F(t))=0. 

Thus, the coupling with the heat bath is described by two 
terms: an operator- valued random force F(t) with mean zero, 
and a mean force characterized by a memory function /j,(t). 
Explicitly, 



n(t) = mjUJj cos(uijt)9(t), 



[6] 



f°° 7re 2 (fcT) 2 na 

U ^ = 3 la ^ W ^ = 3nmc 3 = 3^? ^ ' [3] with the Heaviside step function. Also 



This theoretical result appears to agree with experiment. How- 
ever, as we have previously indicated [12], there are flaws with 
this analysis since: 

(a) missing from (3) is the equipartition term kT/2, which 
is the leading and dominant term, 

(b) radiation damping. The key point is that an atom inter- 
acting with BBR at temperature T is a thermodynamic system. 
Equilibrium is preserved by virtue of the fact that the BBR not 
only gives energy to the atom but also receives energy from 
the atom because of dissipative effects [17, 18], which is a 
beautiful example of the fluctuation - dissipation at work (and 
analogous to Langevin's treatment of Brownian motion, when 
he used the reverse process to introduce a fluctuation force to 
counteract the dissipative force exerted by the fluid). The con- 
clusion is that thermodynamic principles must be used in this 
atomic problem. In particular, the work done in an isothermal 
transition (in this case, the energy supplied by a photon driv- 
ing a transition from the ground state to an excited state) is 
equal to the change in free energy. Thus, it is our basic con- 
tention that the Hollberg - Hall experiment is actually measur- 
ing changes in free energy, as distinct from changes in energy. 
Thus, we next turn to how such changes are calculated. Our 
starting point is the Hamiltonian of the IO system [11] 



F(i)=Y,m^(i), 



[7] 



where q h {t) denotes the general solution of the homogeneous 
equation for the heat-bath oscillators (corresponding to no in- 
teraction). An exact solution can be obtained in the case of an 
oscillator potential V{x) = \Kx 2 = ^moj 2 x 2 , which is best 
displayed as 



x(u) = a(co){F(co) + f(w)}. 



[8] 



Here, the superposed tilde is used to denote the Fourier trans- 
form and a(z) is the generalized susceptibility (response func- 
tion), which is given by 



a(z) = 



izfl(z) 



K 



[9] 



As already remarked, the BBR Hamiltonian is a special case of 
the IO model, for which [11] 



H 



2m 



^7 1 2 / \ 2 

~ m jUj [Qj - x ) 



cf(t). 



[4] 



Here m is the mass of the quantum particle while mj and uij 
refer to the mass and frequency of heat-bath oscillator j. In 
addition, x and p are the coordinate and momentum operators 
for the quantum particle and qj and pj are the corresponding 



Re\fi{uo + i0 + ) 



2e 2 uj< 



[10] 



where the quantity fk is the electron form factor (Fourier trans- 
form of the electron charge distribution). In other words, we 
have allowed the electron to have structure. 

The physically significant results for this model should not 
depend upon details of the electron form factor, subject, of 
course, to the condition that is be unity up to some large fre- 
quency ft and falls to zero thereafter. A convenient form which 
satisfies this condition is 
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fk 



ft 2 



lu 2 + ft 2 ' 



[11] 



Using this in (10), the Stieltjes inversion formula gives 



3c 3 z + ift ' 



[12] 



In addition, we found that the fluctuation force is eE, where E 
is the electric field operator for the free BBR field. The above 
derivation started with the IO model and derived results for 
BBR as a special case. Actually, in our first paper on this sub- 
ject [11], we dealt directly with the BBR Hamiltonian. 

We emphasize again that our QLE given in (5) can be ap- 
plied to many different heat baths of interest (the Ohmic, the 
single relaxation time model, the BBR and so on) but here we 
concentrate on the BBR. It turns out that the BBR model is 
unique in the sense that, as is well known, an essential aspect 
of QED theory is the necessity for mass renormalization. Thus 
the m occurring in the QLE is actually the bare mass and the 
renormalized (observed) mass M is given in terms of the bare 
mass m by the relation [11] 



, , m + 2e 2 n ... 
M = — = m + r e ftM, 



3c 3 



[13] 



where 



AU = -AF , 



[18] 



which is the negative of the result of what we regard as the 
flawed calculation given in (3) above. We conclude that our 
result given in (16) agrees with the results of the Hollberg- 
Hall experiment and that the experiment actually measures free 
energy level shifts. 

3. Equation of motion of a radiating electron 

In a certain sense, the QLE given in (5) [in conjunction with 
(10) and the knowledge that F(t) = eE, where E is the elec- 
tive field operator for the free BBR field] is the required equa- 
tion of motion. However, as already noted, for BBR, the m 
appearing in (5) is the rest mass and thus we must use (13) to 
get the corresponding result in terms of the observed mass M. 
This leads to the result [18] 

(m/fi)z(t) + Mx{t) + V: ff (x) = F eS (t) + f eS (t), [19] 
where 

f eff (t) = f(t) + n-\f(t), [20] 

and similarly for the other "effective" quantities. We note that 
(19) is an exact quantum mechanical result. In the classical 
limit and with V(x) = 0, we obtain 



2e 2 
3Mc 3 



6 x 10~ 24 s. 



[14] 



In the next section, we will return to these results in order to 
express the QLE in terms of M and hence obtain the equation 
of motion of a radiating electron. For now, we point out the 
importance of a(co) in that it leads us to a simple formula for 
F (T), the free energy of the oscillator coupled to the radiation 
field, in the form (11) 



M(ft _1 - T e )'£(*) + Mx{t) = f(t) + ft _1 /(i)- 



[21] 



We note the generality of this result in that we have not yet 
specified the cutoff frequency ft, which, of course, determines 
the form-factor. Also, the principle of causality (response due 
to an external force cannot precede the force) implies that the 
poles of the response function a(u>) must lie in the lower half 
of the complex plane (noting that Imw > [9, 11]) which, in 
turn, implies m > 0. This leads to the conclusion [18] that 



Fo(T} .if mm {^^l}, [15] 

where f(uj,T) is the free energy of a single oscillator of fre- 
quency to, given by 

f(u,T) = fcTlog[l -exp(-hu/kT)] . [16] 

This then led us to the conclusion [11] that the corresponding 
free energy level shift is given by 



ft < r" 1 = 1.60 x 10 23 s~\ [22] 

which rules out the possibility of a point electron [18] and ex- 
plains why the Abraham-Lorentz equation is not an acceptable 
equation for the radiating electron [19]. In essence, we are now 
left with a family of solutions depending on the choice of ft or, 
concomitantly, the choice of electron structure. The simplest 
solution emerges if we choose m = which to equivalent 
to choosing for ft its largest permissible value of r" 1 , corre- 
sponding to choosing the closest approach to a point electron 
consistent with causality. In that case, we obtain 



AF °~ 9M C 2 (fcT)2 ' [17] Mx(t)=f(t)+rJ(t). [23] 

which in 3 dimensions is to be multiplied by 3. It follows from This is rather striking result in that it is only a second-order 

thermodynamics [11] that the corresponding energy level shift equation, it is correct to first order in r e and it is independent 

is of the cutoff frequency ft. We note that the right-side of (23) 
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depends only on the specified external force and thus it is a 
simple equation to solve. Finally, we note that other physically 
reasons choices for the form factor, as distinct from the choice 
given in (12), simply lead to additional higher-order terms in 
(23), such as r^f(t), as is shown explicitly in [18] and [19]. 
Also, we remark that the relativistic generalization of (23) has 
been obtained [20]. 



4. Other miscellaneous applications 

The problem of Brownian motion is a special case of our 
general QLE given in (5) and corresponds to taking f(t) = 0, 
V(x) = 0, and fj,(t - t') = 2m^5(t - t'), so that fl(uj) = m-/ 
where 7 is a constant. In addition, the classical high temper- 
ature limit is also assumed. The end result is that one obtains 
so-called normal (Einstein) diffusion with a diffusion constant 
(kT/mj), which is necessarily a classical result. However, for 
other choices of fl(uj), one obtains anomalous diffusion, of in- 
terest in a variety of applications [21]. In the QED case, for 
which the choice for /i(w)given in (12) is relevant, interest- 
ing quantum effects are manifest. In particular, we find that, at 
T = 0, the result for the diffusion constant contains not only 
quantum effects but also the bare mass appears in the result. In 
addition, it was possible to calculate the spreading of a wave 
packet in a BBR environment [21], a result not amenable to 
calculation by conventional QED methods. 

Next, we point out that our result for the free energy, Fo, 
given in (15), provides the basis for calculating quantum ef- 
fects on the laws of thermodynamics. In that context, we pointed 
out the flaws in a variety of papers which claimed that quantum 
effects could lead to violations of the second and third law of 
thermodynamics [22, 23, 24]. Moreover, we calculated explic- 
itly quantum corrections for various thermodynamic quantities 
(free energy, energy, entropy, and specific heat) for a variety of 
heat bath models [25]. 



5. Conclusions 

Treating the BBR field as a heat bath enabled us to treat 
it in the general context of stochastic physics, with its atten- 
dant powerful results (such as the fluctuation-dissipation theo- 
rem). As a result, we were able to calculate interesting physical 
phenomena not amenable to solution by the conventional tech- 
niques of atomic physics and QED. 
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